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ABSTRACT 

We solve the Riemann problem for the deceleration of an arbitrarily magne- 
tized relativistic flow injected into a static unmagnetized medium in one dimen- 
sion. We find that for the same initial Lorentz factor, the reverse shock becomes 
progressively weaker with increasing magnetization a (the Poynting-to-kinetic 
energy flux ratio), and the shock becomes a rarefaction wave when a exceeds a 
critical value, Uc, defined by the balance between the magnetic pressure in the 
flow and the thermal pressure in the forward shock. In the rarefaction wave 
regime, we find that the rarefied region is accelerated to a Lorentz factor that is 
significantly larger than the initial value. This acceleration mechanism is due to 
the strong magnetic pressure in the flow. We discuss the implications of these 
results for models of gamma-ray bursts and active galactic nuclei. 

Subject headings: active galactic nuclei; gamma-rays: bursts - numerical - MHD 
- relativity 



1. Introduction 

Relativistic jets are believed to exist in active galactic nuclei (AGNs), black hole binaries, 
and gamma-ray bursts (GRBs), but their composition is still poorly understood. It has been 
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argued that magnetic fields could play an important dynamic role in these jets (e.g. Lovelace 
1976; Blandford 1976; Blandford & Znajck 1977; Blandford & Payne 1982; Usov 1992; 
Thompson 1994; Meszaros & Rees 1997; Lyutikov & Blandford 2003; Vlahakis & Konigl 
2003, 2004), but the degree of magnetization, quantified by the magnetization parameter a 
(the ratio of electromagnetic to kinetic energy flux), is poorly constrained by observations. 
GRB afterglow modeling indicates that the ejecta are more magnetized than the ambient 
medium, suggesting a possibly important dynamic role for magnetic fields in GRB jets (Fan 
et al. 2002; Zhang, Kobayashi & Meszaros 2003; Kumar & Panaitescu 2003; Gomboc et al. 
2008). 

A useful diagnostic for the degree of jet-magnetization can be obtained from the inter- 
action between the decelerating jet and the ambient medium. Added magnetic field pressure 
in the jet alters the condition for formation of a reverse shock (RS) as well as its strength 
(Kennel & Coroniti 1984). Analytical studies of the deceleration of a GRB fireball with 
arbitrary magnetization (Zhang & Kobayashi 2005, hereafter ZK05; see also Fan et al. 2004 
for 0" < 1) suggest novel behavior that does not exist in pure hydrodynamic (HD) (cr = 0) 
models (Sari & Piran 1995; Kobayashi et al. 1999). However, consensus on the conditions 
required for the existence of the RS or how Poynting flux is transferred to kinetic flux in the 
interaction region has not yet been achieved (ZK05; Lyutikov 2006; Giannios et al. 2008). 
We present a one-dimensional (1-D) study of the interaction between a magnetized relativis- 
tic flow and a static, unmagnetized external medium. A Riemann problem is solved both 
analytically and numerically over a broad range of a. 



2. The Riemann Problem 

We consider a Riemann problem consisting of two uniform initial states (left and right) 
with discontinuous hydrodynamic properties specified by the rest-mass density p, gas pres- 
sure p, specific internal energy specific enthalpy /i = 1 -|- uj pc? -\- p/pc^, and normal 
velocity . The right state (the medium external to the jet) is assumed to be a cold fiuid 
with constant density, at rest. Specifically, we select the initial conditions: pu — l.Opo, 
Pr — 10~^poC^; '^R = '^K = 0-0) where po is an arbitrary normalization constant (our simula- 
tions are scale-free) and c is the speed of light. The left state (the propagating relativistic 
flow) is assumed to have a higher density and pressure than the right state, as well as a 
relativistic velocity. Speciflcally, pi = lO'^po, pi = l-OpoC^j and = "^2 = 0.995c (7^ ~ 10). 
The fluid is described by an adiabatic equation of state p o: p^ with P = 4/3. 

To investigate the effects of magnetic flelds, we consider a perpendicular fleld component 
in the jet with By = 31.623, 100.0, 316.23, and 447.21 in units of {inpoc^y^^ measured in 
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the laboratory frame, corresponding to a = B"^ /Ari'-f'^hpc^ ~ B'^ / An'-f'^ pc^ being 0.1, 1.0, 10.0, 
and 20.0, respectively. This field is motivated by the predicted toroidal field domination 
at the deceleration radius for GRB outfiows (e.g.. Spruit et al. 2001; ZK05). Increasing a 
increases the total (kinetic plus magnetic) energy density of the left (jet) state. 

3. Results 

We calculate exact solutions of this problem, using the code of Giacomazzo & Rezzolla 
(2006), in the region 0.8 < x < 1.2 with an initial discontinuity at x = 1.0, where x is in 
arbitrary units. 

3.1. Flow-Medium Interaction 

The exact solutions are presented in Figure 1. The four panels display profiles of the 
gas density, gas (and magnetic) pressure, magnetic field strength B^, and the Lorentz factor 
at time t = O.lqj. Different colors represent different a values: 0.1 (black), 1.0 (red), 2.7 
(yellow), 10.0 (green) and 20.0 (blue). The initial Lorentz factor of the left state (jet) is 
7L = 10. 

For cr = 0.1 (black), the solution shows a right-moving fast shock (FS: forward shock; 
S^), a left-moving fast shock (RS: reverse shock; ^S) relative to the contact discontinuity 
(C). In the laboratory frame, the contact discontinuity and the two shocks move to the 
right. 

For o" = 1.0 (red), the solution shows similar profiles {^SCS^) as for a = 0.1. The FS 
is stronger (due to a higher jump in pressure) and slower (more deceleration relative to the 
frame of the contact discontinuity), while the RS is weaker but faster. These features are 
expected from analytical work (ZK05; Giannios et al. 2008), and agree with 1-D relativistic 
MHD simulations (Mimica et al. 2007, 2008). 

When the magnetization of the fiow exceeds a = 2.7, the shock profiles change drastically 
(the significance of this particular value of a is discussed below). For a = 10.0 (green) and 
a = 20.0 (blue), a prominent left-going rarefaction wave {^R) is observed, instead of a 
left-going shock (see also Romero et al. 2005; Mimica et al. 2007). When the rarefaction 
wave propagates into the jet fiow, density and gas pressure decrease, and the fiow velocity 



Here t is in units of x/c with c = 1. 
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Fig. 1. — Profiles of (a) density, (b) gas pressure {solid lines) and magnetic pressure {dotted 
lines), (c) magnetic field (-B^), and (d) Lorentz factor (7) of a = 0.1 (black), 1.0 (red), 
10.0 (green), 20.0 (blue) cases at time t = 0.16. Other parameters: pi = 100.0, pr = 1.0, 
7i = 10.0. The critical value, ctc — 2.7 case is shown as yellow lines. Close-up forward 
shock regions are inserted. Arrows indicate four physically distinct regions: (1) unshocked 
medium, (2) shocked medium, (3) shocked fiow and (4) unshocked fiow corresponding to the 
cr = 0.1 case. 
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increases. The terminal Lorentz factor of the left (jet) state and the FS region reaches 7 ~ 14 
for cr = 10 and 7 > 16 for a = 20. This magnetic acceleration mechanism stems from the 
magnetic pressure in the fio\\EI. 

This magnetic acceleration mechanism is solely a MHD effect and requires the magnetic 
field to generate a rarefaction wave. This is different from the HD/MHD boost mechanism 
proposed by Aloy & Rezzolla (2006), and further investigated by Mizuno et al. (2008) and 
Aloy & Mimica (2008). The HD/MHD mechanism is a purely relativistic mechanism, which 
invokes a relativistic flow perpendicular to the propagation direction of the rarefaction wave. 
The mechanism discussed here occurs even in the Newtonian case, and acts parallel to the 
propagation direction of the rarefaction wave. In general, the acceleration efficiency is smaller 
than that of the HD/MHD boost mechanism (see §3.3 for more discussions). 



3.2. Conditions for Reverse Shock or Magnetic Acceleration 

The magnetic pressure profiles (dotted lines) in Fig. 1(b) reveal the physical condition 
required for the transition from a reverse shock to a rarefaction wave. It is evident that 
the magnetic pressure increases as a increases. In the reverse shock cases (a = 0.1,1), 
the upstream magnetic pressure is lower than the gas pressure in the forward shock, while 
in the rarefaction wave cases {a = 10,20), the upstream magnetic pressure exceeds the 
gas pressure in the FS. Thus, the balance between the upstream magnetic pressure in the 
unshocked flow region and the FS gas pressure in the shocked medium (ZK05; Romero et 
al. 2005) provides the condition separating the two regimes. This condition can be derived 
analytically (see also ZK05). For the interaction between a relativistic flow and an external 
medium, there exist four physically distinct regions: (1) unshocked medium, (2) shocked 
medium, (3) shocked flow, and (4) unshocked flow. Hereafter, Qi denotes the value of a 
quantity "Q" in region "i". From the relativistic shock jump conditions with F = 4/3, 
one can write U2/ p2C^ = (72 ~ 1) — 72 and P2/P1 = 472 + 3 ~ 472. A constant speed 
across the contact discontinuity requires 72 = 73, and the relation between the gas pressure 
and the internal energy gives p2 = M2/3. Thus, the thermal pressure generated in the FS 
region is p2 = (l/3)(72 — 1)(472 + 3)pic^ ~ (4/3)72Pic^. The pressure balance condition is 
Bl/8n'~fl ~ P2, because in pressure balance there is no reverse shock or rarefaction wave, so 
that region 4 and region 3 are matched as -B4 ~ B^, and 74 ~ 73 = 72. Using the definition 



^We note that Romero et al. (2005) and Mimica et al. (2007) also discovered the rarefaction wave regime 
discussed in this paper, but did not investigate the magnetic acceleration mechanism and its astrophysical 
implications in detail. 
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cr = (T4 = Bl/4:7r'ylp4^c^, one can derive a critical (Xc value 

a, = - 74 - 1 474 + 3 — ~ -74— = -7l— • 1 
3 p4 3 P4 3 PL 

The condition for the existence of a reverse shock is a < ac, which is Eq. (31) of ZK05. 
The condition for a rarefaction wave and magnetic acceleration is a > ac- We adopted 
Pi = Pr = 1-0, Pi = Pl = 10^, and 74 = 7l = 10.0, so that the critical value is ac — 2.7. 
Our calculations indicate that ac marks the transition point where neither a reverse shock 
nor a rarefaction wave is established (yellow lines in Fig. 1). To verify this for a larger 
parameter space, we investigate the cx-dependences of various quantities in detail. Fig. 2(a) 
shows the gas pressure in the region through which the reverse shock/rarefaction wave has 
propagated. Initial Lorentz factors are = 5, 10, and 20, respectively. For all cases, we fix 
the flow density at = 10^ and increase B (hence a). The total initial energy density of 
the flow increases with a. In all cases, the gas pressure decreases with a smoothly without 
a sharp transition from the RS regime (solid lines) to the reverse rarefaction wave regime 
(dotted lines). The critical magnetization parameters are ac — 0.7, 2.7, 10.6 for 7^ = 5, 10, 20, 
respectively, derived from the analytical solution Eq.([l]). We notice that in the RS regime, 
the strength of the shock decreases rapidly with increasing a. The critical magnetization 
parameter ac increases with 7^, so that a RS can exist in the high-a regime if 7/, is sufficiently 
large (see also ZK05). 

Another commonly invoked RS-condition states that the shock speed in the fluid frame 
(region 4) is higher than the speed of the Alfven wave, i.e. 7^3 > 7^, where 7^3 j^j^ = 
(74/7rr,rs + 7rr,rs/74)/2, 7^ = (1 + cr)^/^, and 7rs,rr is the Lorentz factor of the reverse 
shock or reverse rarefaction wave in the laboratory frame, which can be calculated using 
the exact solution of Giacomazzo & Rezzolla (2006). Giannios et al. (2008) claim that this 
condition is different from the pressure balance condition. However, in Figures 2(b) and 
2(c), we present the ratios of gas pressure in the FS region to magnetic pressure in the flow, 
Pfs/pb, and 7RSRR/7yi5 flnd that both ratios reach unity at the same critical value ac- 
This suggests that the two reverse shock conditions have the same physical origin, at least 
for the 1-D model studied here. 



3.3. Terminal Lorentz Factor and Magnetic Acceleration Efficiency 

To better understand the magnetic acceleration mechanism, we plot the Lorentz factor 
as a function of a in Fig. 2(d). For the magnetic acceleration case, this is the terminal 
Lorentz factor after acceleration. Because of the dependence of ac on 7^, a higher a is 
needed to achieve acceleration for a higher 7^,. The terminal Lorentz factor can be estimated 
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Fig. 2. — The cr-dependences of (a) gas pressure in the region through which the reverse 
shock (RS; sohd hues) /rarefaction wave (RR; dotted hues) has propagated; (b) the ratio of 
the gas pressure in the forward shock region to the magnetic pressure in the flow; (c) the 
ratio of the Lorentz factor of the propagating reverse shock/rarefaction wave to the Alfven 
Lorentz factor in the rest frame of the fluid; (d) the maximum Lorentz factor in the shocked 
region, in the exact solution. Different initial Lorentz factors have been calculated: 7l = 5 
(black), 10 (red), and 20 (blue). Crosses are the values of the estimated terminal Lorentz 
factor in the 71, = 20 case according to Eq.(2). A constant flow density is adopted, so that 
the total initial energy density of the flow increases with a. 
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Fig. 3. — The dependences of (a) the terminal Lorentz factor 7^; and (2) the acceleration 
efficiency 'jt/lL on the initial Lorentz factor 7^. Solid lines are for the RR regime and dotted 
lines are for the RS regime. Different inital magnetizations have been calculated: a = 10 
(black), 50 (red), and 100 (blue). Crosses are values estimated with Eq.(2) for cr = 100. 

analytically by requiring that the thermal pressure in the FS region balances the magnetic 
pressure in the region through which the rarefaction wave has propagated. For the terminal 
Lorentz factor 7^, this condition can be expressed (roughly) as B^/Srcy^ = (1/3)(7j — l)(47t + 
3)pic^ ~ (4/3)7^^pic^. From the definition of a with B4 ~ B3, this becomes 



Crosses in Fig. 2(d) denote values of terminal Lorentz factors calculated from Eq.(l2]) for model 
parameters, 74 = 7l = 20, pi = Pr = 1.0, and Pi = pi = 10^, in good agreement with the 
exact solution of the Riemann problem in the reverse rarefaction wave regime. 

To investigate the acceleration efficiency, we present in Fig. 3 the terminal Lorentz 
factor 7(, and its ratio to the initial Lorentz factor {'^t/lL) as a function of the initial flow 
Lorentz factor 7/,. While a flow with a higher initial Lorentz factor reaches a higher terminal 
Lorentz factor, a lower initial Lorentz factor implies a higher acceleration efficiency. From 
Eq.([2]) it follows that 74/74 ~ (3crp4/8pi)^/^74 ^^'^ , in good agreement with the exact solution 
of the Riemann problem in relativistic regime. 




(2) 
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4. Summciry and Discussion 

We solved the 1-D Riemann problem for the deceleration of an arbitrarily magnetized 
relativistic flow in a static, unmagnetized medium. For the same initial Lorentz factor, the 
reverse shock becomes progressively weaker with increasing a, turning into a rarefaction 
wave when a > ac, at which point the magnetic pressure in the flow is balanced by the 
thermal pressure in the forward shock. In the rarefaction wave regime, material in the FS 
region is accelerated due to the strong magnetic pressure in the flow. This magnetic acceler- 
ation mechanism may thus play an important role in the dynamics of strongly magnetized, 
relativistic flows. 

Numerical MHD simulations (e.g. Koide et al. 1999, 2000; Nishikawa et al. 2005; Mizuno 
et al. 2007) are essential to understand magnetized relativistic jets. We performed 1-D special 
relativistic MHD simulations of a relativistic flow propagating in an external medium using 
the RAISHIN code (Mizuno et al. 2006). The simulation results are in good agreement with 
the exact solution (Giacomazzo & Rezzolla 2006), serving as a test of the RAISHIN code. 
As MHD simulations can tackle problems for which an exact solution is not known, we plan 
to utilize RAISHIN to solve more realistic configurations (e.g., relativistic shells with a finite 
width and conical geometry, as often envisaged in the GRB problem). 

The magnetic acceleration mechanism discussed here also applies in the Newtonian MHD 
limit. The transition point from a reverse shock to a rarefaction wave is then also given by 
the pressure balance condition, and the terminal velocity of the flow can be estimated from 
the Newtonian shock jump condition as Vt — Csi{p2/pi — l)A/2/r/[(r + l){p2/pi) + (r — 1)], 
where c^i = (Fpi/pi)^/^ is the sound speed in the upstream medium (see also Hawley et al. 
1984). Defining a = (S|/87r)/(p4^^|/2) in the Newtonian hmit, one can derive the terminal 
velocity of the flow, Vt determined by balance between the magnetic pressure and the pressure 
in the forward shock region, p2 ~ B^/SiTT ~ Bl/Sir = ap^vl/^. For a strong shock, (p2 ^ Pi), 
the terminal velocity can be approximated as Vt — a/ ct^Pa/ Pi)/ (F + l)f4. For Vt = v^, one 
derives ac — {T + l)pi/ p^. For F = 5/3, typical for non-relativistic shocks, this expression is 
consistent with Eq.(l) for 7^ = 1. Although the general physics is the same, the dependence 
of the terminal velocity {vt/c)^t is different in the Newtonian and relativistic case. 

Our results have implications for understanding deceleration of strongly magnetized 
outflows, possibly present in GRBs and AGNs. Exact solutions indicate that the condition 
for the existence of a reverse shock is a < cXc, as suggested by ZK05 (cf. Giannios et al. 
2008). The paucity of bright optical flashes in GRBs (e.g., Roming et al. 2006) may, 
among other interpretions, be attributed to highly magnetized flows. Furthermore, the 
magnetic acceleration mechanism discussed here suggests that a and 7 are not independent 
parameters at the deceleration radius. For high-cr flows, the ejecta would experience magnetic 
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acceleration at small radii, before reaching the coasting regime, so that the coasting Lorentz 
factor (i.e., the "initial" Lorentz factor for the afterglow) is at least the "terminal" Lorentz 
factor defined by Eq.([2]). As a result, the high-a and low-7 part of parameter space is 
suppressed. This implies that some region in the ^ — cr parameter spac^ of GRB models 
(ZK05; Giannios et al. 2008) is suppressed as well. Here we only focus on 1-D models with 
Cartesian geometry. Implications for GRB models will be discussed in more detail when this 
Riemann problem is solved in conical jet geometry. 

Variable emission (down to minute timescales) observed in some TeV blazars suggests 
very high Lorentz factors in AGN jets (Aharonian et al. 2007). Models for the production 
of TeV emission appeal to high Lorentz factor jet cores surrounded by lower Lorentz factor 
sheaths (e.g., Ghisellini et al. 2005) or rapid jet deceleration (Georganopoulos et al. 2005) in 
order to reconcile the required jet Lorentz factors with lower Lorentz factors suggested by 
proper motion studies. Our results suggests the possibility of magnetic acceleration occuring 
where highly magnetized jet material overtakes more weakly magnetized jet material. In this 
case the magnetically accelerated Lorentz factor behind the forward shock can significantly 
exceed the Lorentz factor of the overtaken jet material. 

In this study we considerd a static, unmagnetized medium, but for a weakly magnetized 
medium, our conclusions hold as well. However, the conclusions will not apply in the very 
large a regime when the underlying MHD approximation breaks down, e.g., a ~ several 100s 
in the problem of GRBs (Spruit et al. 2001; Zhang & Meszaros 2002). 
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is defined as ■\/ Tdec/fs (Sari & Piran 1995), where rdec is the deceleration radius where the ejecta 
accumulate from the external medium a mass 7q"^ times their own mass, and r^. is the spreading radius 
where the width of the ejected shell starts to increase due to propagation of a sound wave. 
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